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We report results from an updated study of the suppressed decays B~ DK~ and B~ — > D*K~ 
followed hy D ^ K+iv', where D<^*^ indicates a D<'*^° or a 5^*'° meson, and D* Div° or 
D* Df. These decays are sensitive to the CKM unitarity triangle angle 7 due to interference 
between the 6 — ^ c transition B~ D^*^°K~ followed by the doubly Cabibbo-suppressed decay 

7?" — > K'^n^, and the b ^ u transition B^ — > D^*^^K^ followed by the Cabibbo-favorcd decay 
— >■ K^%^ . We also report an analysis of the decay B^ — >■ D^'V" with the D decaying into the 
doubly Cabibbo-suppressed mode D — ^ K'^t^^ . Our results are based on 467 million T(AS) — >■ BB 
decays collected with the BABAR detector at SLAC. Wc measure the ratios T?.*-*-* of the suppressed 
{[K'^-k~]dK~ /-k') to favored {\K^-k^]dK~ /-k') branching fractions as well as the CP asymmetries 
A^*'' of those modes. Wc sec indications of signals for the B^ — > DK~ and B^ — > ^*dttO^~ 
suppressed modes, with statistical significances of 2.1 and 2.2a, respectively, and we measure: 

Udk = (1.1 ± 0.6 ± 0.2) X 10"^ Adk = -0.86 ± 0.47 t^le, 

■^(i37r0)jf = (1-8 ± 0.9 ± 0.4) X 10"^ AId^o)k = +0-77 ± 0.35 ± 0.12, 

'R-*D-y)K = (1-3 ± 1-4 ± 0.8) X 10-^ A*D-y)K = +0.36 ± 0.94 tUl, 

where the first uncertainty is statistical and the second is systematic. We use a frequentist approach 
to obtain the magnitude of the ratio rs = \A{B- D°K-)/A{B- D°K-)\ = (9.5l|;J)%, 
with rs < 16.7% at 90% confidence level. In the case of B' D* K' we find r|j = \A{B- 
D*°K-)/A{B- D*°K-)\ = (9.6ig;5)%, with < 15.0% at 90% confidence level. 
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PACS numbers: 13.25.Hw, 14.40.Nd, 12.15.Hh, ll.30.Er 



I. INTRODUCTION 

The Standard Model accommodates CP violation 
through a single phase in the Cabibbo-Kobayashi- 
Maskawa (CKM) quark mixing matrix V [ijj. In 
the Wolfenstein parameterization 0, the angle 7 = 
arg {—VudV*ij/VcdV*,j) of the unitarity triangle is related 
to the complex phase of the CKM matrix element Vub 
through Vub = iKible"*'''. A theoretically clean source of 
information on the angle 7 is provided by B~ — D^*^ 
decays, where D^*"^ represents an admixture of D^*'>'^ and 
D'-*^^ states. These decays exploit the interference be- 
tween B- D(*'>°K- and B' -^_&*^°K- (Fig.[T]) that 
occurs when the and the decay to common 

final states. 



where r r = \A{B- D°K-)/A{B- D°K-)\ w 
10% [Jd and ro = \A{D^ K+tt-)/A{D° 
K-TT+)\ = (5.78 ± 0.08)% [8] are the suppressed to fa- 
vored B and D amplitude ratios. The rates in Eq. ([1]) 
depend on the relative weak phase 7 and the relative 
strong phase S = 5b + between the interfering am- 
plitudes, where 5b and 5d are the strong phase differ- 
ences between the two B and D decay amplitudes, re- 
spectively. The value of 5d has been measured to be 
5d — (201.9122.4)° M, where we have accounted for a 
phase shift of 180° in the definition of 5b between Ref. Q 
and this analysis. 

The main experimental observables are the charge- 
averaged decay rate and the direct CP asymmetry, which 
can be written as 




\J 



£)(*)() 




T^DK 
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DK 



< D 



2 tbtd cos 7 cos ^ 



(2) 



FIG. 1: Feynman diagrams for B~ d'^')o x'-*''- and 
j)(*)o^(*)- rpj^g latter is CKM and color-suppressed with 
respect to the former. 



In the Atwood-Dunietz-Soni (ADS) method [3], the 
D° from the favored c amplitude is reconstructed 
in the doubly Cabibbo-suppressed decay K~^Tr~ , while 
the D° from the u suppressed amplitude is recon- 
structed in the favored decay K~^Tr~. The product 
branching fractions for these final states, which we de- 
note as [K+7:-]dK^ (P- DR-) and [K+t:-]d'K- 
{B~ D*K^), are small (~ 10"''), but the two inter- 
fering amplitudes are of the same order of magnitude, 
and large CP asymmetries are therefore possible. The 
favored decay mode B~ — )■ [K~tt^]bi{,)K~ is used to 
normalize the measurement and cancel many systematic 
uncertainties. Thus, ignoring possible small effects due 
to D mixing and assuming no CP violation in the nor- 
malization modes, we define the charge-specific ratios for 
_B+ and B^ decay rates to the ADS final states as 



n 



DK 



2 tbTd cos(±7 + 5), 



(1) 



Adk 



n 



7^ 



DK ^ '^DK 

2 TBro sin J sm5/7lDK- (3) 



The treatment for the D*K mode is identical to the DK 
one, but the parameters and 5b are not expected to be 
numerically the same as those of the DK mode. Taking 
into account the effective strong phase difference of tt 
between the D* decays to D7 and -Dtt^Q, we define the 
charge-specific ratios for D* as: 



n 



*± 

(Dtt«)K 



_ T{[K^t:^]d._,d^oK^) 

- T{[K±n^]D'^D.oK±) 

= r*i +rl + 2r*BrDCOs{±-i + 5*), (4) 



7^ 



(D-i)K 



_ T{[K^n^]D.^D^K^) 
T{[K±n^D'^D-,K±) 



-2r|jri5COs(±7 + <5*), (5) 



with r*g 



\A{B~ ^ D*"K-)/A{B- D*°K~)\ and 
i5* = 5%+5d, where 5'^ is the strong phase difference be- 
tween the two B decay amplitudes. The charge averaged 
ratios for D* — !• Dtt'^ and D* Dj are then: 
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Definitions of tlie direct CP asymmetries -4*^^o-|;^ and 
follow Eq. ©. 
This paper is an update of our previous ADS analy- 
sis in Ref. [1], which used 232 x 10^ BB pairs and set 
90% C.L. upper limits Hdk < 0.029, 7^*^^o)^ < 0.023 
and Ti-*jj^^i( < 0.045. In addition to an increased data 
sample, new features in the analysis include a multi- 
dimensional fit involving the neural network output used 
to discriminate the signal from the continuum back- 
ground, rather than a simple cut on this variable as was 
done in the previous analysis. We also include measure- 
ments of the ratios of the doubly Cabibbo-suppressed to 
Cabibbo- favored Z3^*V decay rates, 

and of the corresponding asymmetries. These measure- 
ments are used as a check for the B~ — [K^-k~]jj{,)K~ 
ADS analysis. In the D'^*^'k case, we expect that the ratio 
_^{*){D-n) ^-^^ amplitudes is suppressed by a 

factor |T4(iKis/^tidVcs| compared to the D'^*^K case, if we 
assume the same color suppression factor for both decays. 
One expects therefore r^''^'^'^^ r^'' x t an'^ 6'c ~ 5 x 
10"'^ '^ro, where 9c is the Cabibbo angle and where we 
have assumed r^-* = 10%. Neglecting higher order terms, 
n^*l ~ r% and A^^l - 2r^*Han2 6*^ sin7 sin(5(*)/rD. 
Hence, the maximum asymmetry possible for D^*^t: ADS 
decays is 2r^^ tan^ Oc/td ~ 18%. 

II. THE BABAR DETECTOR AND DATASET 

The results presented in this paper are based on 
467 x 10^ T{AS) BB decays, corresponding to an inte- 
grated luminosity of 426 fb^^ (on-peak data). The data 
were collected between 1999 and 2007 with the BaBAR 
detector [l^ at the PEP-II e+e~ colhder at SLAG. In 
addition, a 44 fb~^ data sample, with center-of-mass 
(CM) energy 40 MeV below the T{AS) resonance (off- 
peak data), is used to study backgrounds from continuum 
events, e"'"e~ qq {q — u, d, s, or c). 

The EaBAR detector response to various physics pro- 
cesses as well as to varying beam and environmental con- 
ditions is modeled with simulation software based on the 
Geant4 [11] tool kit. We use EVTGEN [l^ to model the 
kinematics of B meson decays and JETSET [l^ to model 
continuum processes e~^e~ — > qq. 

III. ANALYSIS METHOD 

A. Basic Requirements 

We reconstruct B- D^*^K- and D^*\- 
with the D decaying to K'lr^ (right-sign (RS) de- 
cays) and K^TT~ (wrong-sign (WS) decays). For de- 



cays involving a D*, both D* — )• Dtt" and D* — )• 
modes are reconstructed. Charged kaon and pion candi- 
dates must satisfy identification criteria that are typically 
85% efficient, depending on momentum and polar angle. 
The misidentification rates are at the few percent level. 
We select D candidates with an invariant mass within 
20 MeV/ (about 3 standard deviations) of the known 
mass [1J|. All D candidates are mass and vertex 
constrained. For modes with D* Dvr" or D* D^, 
the mass difference Am between the D* and the D must 
be within 4MeV/c^ (~ 4(t) or 15MeV/c^ (~ 2a), respec- 
tively, of the nominal mass difference [l^ |. 

For the WS decays B^ — >■ [K^Tr^]DK'^, two impor- 
tant sources of background arise: the first from B^ — >■ 
[tt^ K^]]:)K^ (in which the K and tt in the D decay 
are misidentified as tt and K) and the second from 
B^ [K^K^]]jTr^ (when the tt^ pair has an in- 
variant mass within 20MeV/c^ of the nominal mass). 
To eliminate the first background, we recompute the in- 
variant mass (Mswitch) of the h^h'^ pair in h^h'~ 
switching the mass assumptions on the h'^ and the h'~ . 
We veto candidates with Afgwitch within 20 MeV/ of the 
mass [l3|. To eliminate the second background, we 
also veto any candidate where the KK invariant mass is 
within 20 MeV/ c'^ of the D° mass. To ensure the same 
selection efficiencies, these criteria are applied both to 
B± ^ [K^Tr^]oMK^ and to B± ^ [if±7r=F]^(.)if± 
candidates. These veto cuts are 88% efficient on signal 
decays. 

We identify B candidates using two nearly inde- 
pendent kinematic variables that are customarily used 
when reconstructing B-meson decays at the T(4S'). 
These variables are the energy-substituted mass, jties = 

y^(f + Po • PbY /Eq — p^g and energy difference /S.E = 

Eb ~ 5%/^' where E and p are energy and momentum, 
the asterisk denotes the CM frame, the subscripts and 
B refer to the T{AS) and B candidate, respectively, and 
s is the square of the CM energy. For signal events 
Toes = TT^B^ [3 and AB = within the resolutions 
of about 2.6 MeV/c^ and 17 MeV, respectively. We re- 
quire that all candidates have | ABj < 40 MeV and we use 
Toes in the fit to extract the number of signal events. 

The average number of S — >■ D^*'>K candidates recon- 
structed per selected event is about 1.4 in B — >• DK sig- 
nal Monte Carlo (MC) events and about 2 for B -> D*K 
signal MC events. This is mostly due to the cross- 
feed between the DK and the D*K final states. For 
all events with multiple B — D^*^ K candidates, we re- 
tain only one candidate per event, based on the smallest 
value of |A_E|. This method does not bias the sample 
since Ai? is not used to extract the number of signal 
events. After this arbitration, less than 0.4% (0.5%) of 
the B DK {B D*K) signal MC events selected 
are reconstructed as B ^ D*K {B -J> DK). About 10% 
of the B — )• D*jj^oK events selected are reconstructed as 
B D*^^K and about 2% of the B D*jj^K events 
selected are reconstructed as i? D*nnK. 
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TABLE I: Selection efficiencies, after correction for known 
data/MC differences, for — > [K^Tv^]jj(,)h^ i^ws) and 
— > [K^7v^]jj(,)h^ (ehs), and efficiency ratio ews/^RS- 



Cliannel ews (%) ens (%) ews/^Rs (10^ 

DK 26.5±0.1 26.6±0.1 99.6±0.5 

^D^oK 13.3±0.1 13.2±0.1 100.6±1.1 

D*DyK 17.4±0.1 17.5±0.1 99.8±0.8 

Dn 26.0±0.1 26.5±0.1 97.9±0.5 

D'jy^oTv 14.3±0.1 14.8±0.1 96.4±0.9 

D'n^n 18.8±0.1 19.5±0.1 96.3±0.7 



The B — 7> D^*^Tr analysis is performed independently of 
the B — > D^*^K analysis, but uses the same multiple can- 
didate selection algorithm. A summary of the selection 
efficiencies for the WS modes [K^n^]DMh^ {h=K,TT) 
and the RS modes [Tr^iiT^J^x.) /i^ is given in Table HI 

B. Neural Network 

After these initial requirements, backgrounds domi- 
nantly arise from continuum events, especially e^e^ — > 
cc, with c ^ D°X, W ^ K+TT- and c -> D°X, 
— >■ if'+anything. The continuum background is re- 
duced by using neural network techniques. To select the 
discriminating variables used in the neural network, we 
rely on a study performed for the previous version of this 
analysis [1], and we consider the seven quantities listed 
below: 

1. Two event shape moments Lq — J2iPi^ L2 = 

Pi cos^ 9i, calculated in the CM frame. Here, pi 
is the momentum and 9i is the angle with respect to 
the thrust axis of the B candidate; the index i runs 
over all tracks and clusters not used to reconstruct 
the B meson (rest of the event). These variables 
are sensitive to the shape of the event, separating 
jet-like continuum events from more spherical BB 
events. 

2. The absolute value of the cosine of the angle in 
the CM frame between the thrust axes of the B 
candidate and the detected remainder of the event, 
I cos 9t I • The distribution of | cos 9t \ is approxi- 
mately uniform for signal and strongly peaked at 
one for continuum background. 

3. The absolute value of the cosine of the CM angle 
between the B candidate momentum and the beam 
axis, I cos 1- In this variable, the signal follows a 
1 — cos^ 9b distribution, while the background is 
approximately uniform. 

4. The charge difference AQ between the sum of the 
charges of tracks in the D*-*-' hemisphere and the 
sum of the charges of the tracks in the opposite 



hemisphere, excluding the tracks used in the recon- 
structed B, and where the partitioning of the event 
into two hemispheres is done in the CM frame. 
This variable exploits the correlation occurring in 
cc events between the charge of the c (or c) in a 
given hemisphere and the sum of the charges of all 
particles in that hemisphere. For signal events, the 
average charge difference is (AQ) = 0, whereas for 
the cc background (AQ) « | x Q^, where Qs is 
the charge of the B candidate. 

5. The product Qb ■ Qk, where Qk is the sum of the 
charges of all kaons in the rest of the event. In 
many signal events, there is a charged kaon among 
the decay products of the other B in the event. 
The charge of this kaon tends to be highly corre- 
lated with the charge of the B. Thus, signal events 
tend to have Qb ■ Qk < — 1- On the other hand, 
most continuum events have no kaons outside of 
the reconstructed B, and therefore Qk ~ 0. 

6. A quantity M.KI, defined to be zero if there are 
no leptons (e or /i) in the event, and, if a lep- 
ton is found, taken to be equal to the invariant 
mass of this lepton and the kaon from B (bache- 
lor K). This quantity differentiates between con- 
tinuum background and signal because continuum 
events have fewer leptons than BB events. Further- 
more, a large fraction of leptons in cc background 
events are from D — >■ Kiv^ where the kaon becomes 
the bachelor kaon candidate, so that the average 
Mki in cc events is lower than in B signal events. 

7. The absolute value of the measured proper time in- 
terval between the two B decays, |Ai|. This is cal- 
culated from the measured separation, Az, between 
the decay points of the reconstructed B and the 
other B along the beam direction, and the known 
Lorentz boost of the initial e+e~ state. For con- 
tinuum background, \/S.t\ is peaked at 0, with most 
events having \/S.t\ < 2ps, while it is less peaked 
and can extend beyond 5ps for B"^ — )• D^*^h'^ sig- 
nal events. 

The neural network is trained with simulated contin- 
uum and signal [K^Tr^]j^(,)K^ events. Only wrong-sign 
D^*'>K candidates are used in the training, but the neural 
network is used in the analysis of all the D'^*^K^ channels. 
The distributions of the neural network output (NN) for 
signal-enriched right-sign control samples are compared 
with expectations from the MC simulation in Fig. ^a.) 
(DK) and Fig.[2rd) (Dtt). The agreement is satisfactory. 
In the same figure, the NN spectra of background control 
samples (off-peak data) are compared with expectations 
from continuum qq MC. Since we do not expect these dis- 
tributions to be exactly the same for the right-sign and 
wrong-sign background samples, they are shown sepa- 
ratelyforthe [K^TT^]j,^,,K^{Fig.m.b)), [K^n^]o,,,K^ 
(Fig.Hc)), [K^7T^]aM7r^{Fig.m.e)) and [ifT^±]^,.,^T 
(Fig. mf)) channels. To increase the statistics, the toes 
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FIG. 2: (color online). Signal and background distributions of the neural network output, and results of the NN verifications 
for DK (a), D^*^ K (b, c), Div (d) and D'*'7r (e, f) candidates. (a,d): Dh^ right-sign candidates, signal-enriched by a cut on 
the AE, mES signal region. Shaded plain histograms are MC expectations for qq background (dark gray/blue), bb background 
(middle gray/green) and — >■ Dh^ signal events (light gray/yellow). Points with error bars are on-peak data. (b,e): D^*''h^ 
wrong-sign background. (c,f): D^*^h^ right-sign background. Plots b, c, e, and f are normalized to unity. The dotted line 
histograms show the distribution of simulated continuum events. The off-peak data used to check the A'^A'^ are overlaid as data 
points. To increase the statistics, the tties and AE requirements on the off-peak and continuum MC events have been relaxed, 
and Dh^ and D*h^ contributions have been summed. 



and Ai? requirements on the off-resonance and contin- 
uum MC events have been relaxed, and the Dh^ and 
D*h^ contributions have been summed, after checking 
that they are in agreement with each other. Good agree- 
ment between data and the simulation is observed in 
all channels. Good agreement between the D'^*^ K and 
the D^*'>'K background NN distributions is also visible in 
Fig. [21 while on the contrary the background NN distri- 
bution of wrong-sign decays is clearly different from the 
background NN distribution of right-sign decays. We 
have examined the distributions of all variables used in 
the neural network, and found good agreement between 
the simulation and the data control samples. Finally, 
we examined the TV TV distributions in the signal MC for 
the different B signal channels, right-sign and wrong-sign 
separately {Dn, D*7t, DK, D*K) and did not observe 
any significant difference between these channels. 



C. Fitting for event yields and 7?.'*' 

The ratios Ti-^*^ are extracted by performing extended 
unbinned maximum likelihood fits to the set of variables 
Toes, NN, and /sign, where /sign is a discrete variable 
equal to for WS events and to 1 for RS events. We 
write the extended likelihood C as 



N 



L = 



N\ 



where the vector x indicates the variables (toes, NN , and 
/sign) and 9 indicates the set of parameters which are 
fitted from the data. N is the total number of signal and 
background events, and N' = Ni is the expectation 
value for the total number of events. The sum runs over 
the different signal and background categories i which 
will be detailed below. The probability density function 
(PDF) /(xj|6') is written as the sum over the different 
signal and background categories 



/(x,|0,A^') 



N' 



where /i(x|6') is the product F(toes) x G{NN) x /f(/sign) 
of an Toes component /"(toes), a NN component G{NN) 
and a two-bin histogram //(/sign) set to (1,0) for the WS 
category and (0,1) for the RS category. The NN dis- 
tributions are all modeled by histograms with 102 bins 
between -1.02 and 1.02. 

The fits are performed separately to each of the Dtt, 
D*D^oTT, D}y^n, DK, D*jy^oK and D}j^K samples. They 
are configured in such a way that TZ^*^ is an explicit fit pa- 
rameter: for the B signal, we fit for the number of right- 
sign decays N^s and the ratio TZ'-*^ = Nws/{c x Nns), 
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where Nws is the number of wrong-sign signal events 
and c is the ratio of the wrong-sign to right-sign selection 
efficiencies. For B D^*^K, the factor c is consistent 
with unity within the statistical precision of the simu- 
lation (Table |l| and is set to this value in the fits. For 
B — > D^*'>7r, c differs slightly from unity due to differ- 
ent particle identification cuts applied at an early stage 
of the event selection and we use therefore the values of 
Table HI in the fits. 

The following signal and background categories are 
used to describe each sample in the fits: 

1. The right-sign signal B^ [K^tt^]jj(,)K^ /tt^: 
its TOES spectrum is modeled by a Gaussian func- 
tion Gsigi'tn-Es) whose mean and width are deter- 
mined from the fit to data. The NN PDF AfAfsig 
is constructed from the NN spectrum of the B^ ^■ 
Dh- signal MC. 

2. The wrong-sign signal B^ [K^tt^]i,(,)K^ /n^: 
its Toes and NN spectra have the same 
parametrizations Gsig {'tnEs) and AfAfgig as the 
right-sign signal. 



3. The right-sign combinatorial background from qq 
{q — u,d,s,c) events into [K~'!t~^]K~ {DK) or 
[if "Tr+jTT^ (Dtt): its toes component is modeled 
with the ARGUS function [l^ Aqq{mEs) whose 
shape and endpoint parameters, (^qq and toq, are al- 
lowed to vary in the fit. The NN PDF A/'A/'^^"^^ is 
constructed from the NN spectrum of [K~7r^]K^ 
{DK) or [Ar^Tr+Jvr^ (Dtt) candidates in the qq con- 
tinuum MC (Figs. [2t and [2f), where the AE re- 
quirement has been extended to \AE\ < 200 MeV 
and the DK and D*K (or Dtt and D*Tr) samples 
have been summed to increase the statistics. 

4. The wrong-sign combinatorial background from qq 
events into [K+Tr^]K^ {DK) or [K+TT^]n^ {Dtt): 
its Toes component is parameterized by the same 
ARGUS function Aqq{mEs) used for the right-sign 
component. The A^ A^ PDF MM^^^^ is constructed 
from the A^Af spectrum of [K+t:-]K- {DK) or 
[A'+7r~]7r~ (-Dtt) candidates in the qq continuum 
MC (Figs. [3) and Wl- 

5. The right-sign combinatorial background from BB 
events into [K~tt+]K^ {DK) or [if^vr+jTr" {Dtt)-, 
excluding the peaking background which is consid- 
ered in category 7: its toes component is described 
by an ARGUS function [l^ A^g^^ (toes) with shape 
parameter Cs^^"* fixed to its value determined from 
BB MC, after removal of the B -J> D^*) AT/tt signal 
events. The A^A^ PDF used to describe this back- 
ground is the PDF AfAfsig describing the A^A^ spec- 
trum_of the B~ D^*^h~ signal MC. The number 
of BB right-sign combinatorial background events 
is allowed to vary in the Dh~ fits but is fixed to 
the MC prediction in the D*h~ fits (see below). 



6. The wrong-sign combinatorial background from 
BB events into [K+n~]K- {DK) or [AT+Tr-jTr" 
{Dtt), excluding the peaking background which is 
considered in category 8: its toes component is de- 



scribed by an ARGUS function [15| ^^"^^(toes) 



with shape parameter Cs^^'' fixed to its value de- 
termined from the BB MC, after removal of the 
B D(*)A'/7r signal events. The A^A^ PDF used 
to describe this background is the PDF NN sig de- 
scribing the A^A spectrum of the B^ D^*'>h^ 
signal MC. The number of BB wrong-sign combi- 
natorial background events is allowed to vary in the 
Dh'' fits but is fixed in the D*h^ fits (see below). 



7. The background from BB events in the right-sign 
component peaking in toes inside the signal re- 
gion (peaking background) : this background is dis- 
cussed in more detail in Section ITVl For the DK^, 
Dtt^ and ZJ^^qAT^ categories, the peaking part 
of the BB background toes spectrum is described 
by the same Gaussian function Gsig{'m'Es) as the 
signal. This component is therefore indistinguish- 
able from the signal and its rate has to be fixed 
to the MC predictions. For the D^^qTt^, D*j-i^tt^ 
and the D*^^K'^ categories, the TOes component is 
described by an asymmetric Gaussian whose shape 
parameters and amplitude for each category are de- 
termined from a fit to the toes spectrum of BB MC 
events, after vetoing the B^ — )• D'^*^h^ signal com- 
ponent. For all categories, the AA PDF used to 
describe this background is the PDF NN sig de- 
scribing the A^A^ spectra of the B D^*^h^ signal 
MC. 



8. The peaking background from BB events in the 
wrong-sign component: the treatment is similar 
to the previous component but G sig {'mEs) is used 
to describe the toes spectrum of the DK"^, Dtt"^, 
D*jy^oK'^ and D*j-,^K'^ categories, while an asym- 
metric Gaussian is used to describe the toes spec- 
trum of the D*jj^oTt'^ and DJ^^tt^ categories. 

To summarize, we fit for the number of right-sign sig- 
nal events A^as, the ratio TZ = Nws/{c x Nrs) of wrong- 
sign to right-sign events, the number of wrong-sign and 



and 



right-sign qq combinatorial background events, N^'^g 
N^P , and for Dh^ the number of wrong-sign and right- 



sign BB combinatorial background events, A^, 



(BB) 
WS 



and 



A]jg ' . We fix to their MC expectations the numbers 
of wrong-sign and right-sign BB peaking background, 

^ws '^^'' ^'^'^ N^^^'^^\ as well as the number of BB 
combinatorial background events for D*h^. The other 
parameters fitted are the reconstructed toes peak and 
resolution, mg and (Jms i and the qq continuum back- 
ground shape parameter and endpoint, (^qq and toq. 
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TABLE II: Charmless background channels and branching 
fractions, Dh"^ channels affected by this background and 
background yields expected in our data sample. 







Affected 




Estimated 


Mode 




channel 


B(io-«) 


Yield 


K- n+ 


TT 


Dn RS 


55 ± 7 lU] 


67.1±9.7 


K+ TT- 


7r~ 


Dn WS 


< 0.9 [16] 


<1.1 


K- 7r+ 


K- 


DK RS 


< 0.2 [16] 


<0.2 


K+ TT" 


K~ 


DK WS 


5.0 ±0.7 \17] 


6.0±0.8 



IV. STUDY OF BB BACKGROUNDS 

We study the BB background for each signal cat- 
egory {Dtt, D*tt DK, D*K) and charge combination 
(right-sign and wrong-sign) using a sample of e+e~ — ?> 
T(45) — ^ BB MC events corresponding to about 3 times 
the data luminosity. In addition, dedicated Monte Carlo 
signal samples are used to estimate the background from 
B~ — > Dh~ events and the background from the charm- 
less decay B~ — K^tt~ K~ . We identify three main 
classes of background events which can peak in toes in- 
side the signal region and mimic the D^*^7r and D'^*'>K 
signal: 

1. Charmless B decays B^ h^h^h^ [h = tt,K): 
we list in Table HIl the 3-body charmless decays af- 
fecting our analysis, their branching fractions [l3] 
and the numbers of reconstructed events expected 
in the affected modes after the selection. Due to 
the particle identification criteria used in the anal- 
ysis only decays with the same final state parti- 
cles as our signal modes contribute significantly to 
the background. These events are indistinguish- 
able from the Dh^ signal if the K^n^ invariant 
mass is consistent with the D mass. The two de- 
cays affected by a significant charmless background 
are right-sign B^ — > [i^^Tr+J^iTr^ and wrong-sign 
B- [K+tt-]dK-. JQsing B' K-tt+tt- 
events selected in the BB Monte Carlo sample, we 
estimate the efficiency of B^ K''tt^-k~ events 
to be reconstructed as a [if^Tr+JuTr" candidate 
as (0.26 ± 0.02)%. The corresponding background 
is estimated to be 67.1±9.7 events, where the er- 
ror is dominated by the statistical uncertainty on 
the B^ K^TT^TT^ branching fraction. The 
efficiency of B^ — >■ K~^Tr~K~ events to be re- 
constructed as [K'^tt~]dK~ WS candidates is de- 
termined from a high statistics dedicated B^ ^■ 
K~^iT~K~ signal Monte Carlo sample, and is found 
to be (0.27 ± 0.01)%. The corresponding peaking 
background from B^ — > K^t:^ K~ events mimick- 
ing B~ — [K^'k—]dK~ WS decays is estimated to 
be 6.0±0.8 events, where the error is dominated by 
the statistical uncertainty on the B~ — >■ K'^tt~K~ 
branching fraction. From a fit to data selected in 



the D mass sidebands, we cross-check this predic- 
tion and find 6.5±4.0 peaking events, in good agree- 
ment with the MC prediction. We also check that, 
because of the tight Am cut applied to the D* de- 
cay products, the B~ — > D*h^ channels are not 
affected by charmless peaking backgrounds. 

2. Events of the type B^ Dh^: this background 
is estimated by running the analysis on a sample 
of B~ — Dh~ signal MC events properly renor- 
malised to the data sample, and fitting the toes 
spectra of the selected events to the sum of a Gaus- 
sian signal and a combinatorial background. We 
find that a peaking background of 2.6 ± 0.4 events 
is predicted in the B~ — >• [K^n—]DK~ WS chan- 
nel. This component is dominated (2 events out of 
2.6) by decays B~ [K-K+]DTr~ failing the D 
mass veto and by WS decays B~ — )• [K^tt~]d't~ 
where the tt" is misidentified as a K~ . For the 
D*K channels, the B^ — ^ [K~ K+]dtt~ contribu- 
tion is suppressed by the Am cut on the D*-D 
mass difference, and the WS D*tt contribution is 
0.5±0.1 events for D* Dn" and 0.6±0.2 events 
for D* — > Dj. Another background of the same 
type occurs in the right-sign DK decays. It consists 
of events B~ — ^ [if ~7r+]^(,)7r~ where the bachelor 
TT~ is misidentified as a K^ , which fake the RS 
signal B~ — > [K~Tr^]jj{t)K~ . This contribution is 
predicted by the simulation and has been verified 
in the data by fitting the AE spectrum of D^*^K 
candidates in the toes signal region, which shows 
a second peak due to candidates, shifted by 
50 MeV with respect to the signal. 

3. Other decays: this component is estimated by fit- 
ting the Toes spectra of BB MC events, after re- 
moving the charmless and B~ — >■ Dh~ components. 
For B^ — > [K^tt^]dK^ WS decays, the peaking 
component is estimated to be 4 ± 3 events, where 
the uncertainty is dominated by the statistical er- 
ror on the simulated data. The main sources of 
peaking background which could be identified are 
listed in Table |llll They include B° D*+h- re- 
constructed as B^ — > D*"h~, semi-leptonic decays 
B° D**-e+Ve {D**- £)(*)07r-, D° K+tt-) 
where the e"*" is missed, faking the WS signal B^ — > 
[_ft'+7r~]£,(,)7r~, and decays B~ D^*^p~ faking 
the RS signal i?^ — >■ [K~Tr^]jj(,)Tr~ . 

A summary of the BB background studies is given in 
Table Uni for B D'^*^n and B -> D'^*^K. For each 
channel, the toes spectra of events selected in the BB 
MC simulation (after removing the corresponding signal) 
were fitted by the sum of a combinatorial background 
component and a peaking component, using the same 
parametrization described in Sec. IIII CI The average 
number of BB combinatorial and peaking background 
events predicted by the simulation are given in Table Hill 
together with the main sources of peaking events and the 
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TABLE III: Expected numbers of signal and BB background events, peaking background parametrization and dominant sources 
of peaking backgrounds for B — > D^*^n and B — > D^*^ K. ]\['^'^^'> jg the combinatorial part of the background, parametrized by 
an ARGUS function, and A'^^^"'"' is the component peaking in rriES , parametrized by either a Gaussian function or a bifurcated 

Gaussian function. The average event yield expected for the WS signal is computed assuming r^'' — 10% and no interference 
term (cos 7 x cos 5 = 0). 



Mode 


Signal yield 


7 Tr\ 

pr(comb) 

BB 


7 n 

BB 


Peaking bkgd. parametrization Peaking bkgd. 


sources 


Dn~ WS 


86 


93.7±6.0 


10.6 ± 3.0 


Gaussian 






Dh.on- WS 


31 


24.7±8.3 


29.0 ±8.7 


Bifurcated Gaussian 






Di,.,TT- WS 


25 


111±9 


47 ±7 


Bifurcated Gaussian 


Dl~e'^Ue, D{' 


e+i/,, and D^*'>°p° 


Dn- RS 


24240 


307.3±11.7 


222.0±10.3 


Gaussian 


K~TV'^TV~, [cc] 


iR- 




8931 


620.7±33.7 


507.3 ± 33.3 


Bifurcated Gaussian 


D*p~, D*+7r" 




D})^Tr~ RS 


7242 


1225±64 


2432 ± 67 


Bifurcated Gaussian 




, and -Dq^oTt" 


DR- WS 


26.3 


107.0±6.3 


12.6 ±3.1 


Gaussian 


Dh-, R-K+Tx- 


D'^^oK- WS 


8.5 


17.3±2.7 


2.7 ± 1.6 


Gaussian 






DhjK- WS 


6.8 


68.3±5.3 


6.0 ±2.4 


Gaussian 






DK' RS 


1944 


50.7±5.3 


299.3±10.7 


Gaussian 






D'^^oK- RS 


618 


56.0±6.7 


127.0 ±8.3 


Gaussian 






Di,^K- RS 


503 


66.0±14.7 


326.7 ± 17.3 


Bifurcated Gaussian 




R- 



functional shapes chosen to describe the peaking back- 
ground. The numbers of signal events expected are also 
given for comparison. For the B — >■ D*K WS channels, 
we could not identify a specific source of peaking back- 
ground due to the lack of statistics in the simulation. For 
all channels, we use the values of the peaking components 
summarized in Table lllll in the maximum likelihood fit. 
Statistical uncertainties in the expected yields are incor- 
porated in the corresponding systematic uncertainties. 

V. RESULTS 
A. Results for B D^'K 

The results for B TT are displayed in Fig. [3] 

(right-sign modes) and Fig. |4] (wrong-sign modes). They 
are summarized in Table IIVI Clear signals are observed 
in the B — ^ Dtt and in the B — > D*jj^oTt WS modes, 
with statistical significances of la and 4.8cr, respectively. 
The significance is defined as -^Z— 21n(£o/'Cmax), where 
'Cmax a-nd Cq are the likelihood values with the nominal 
and with zero WS signal yield, respectively. For B — )■ 
D}j^n WS decays, the significance is only 2a, due to the 
large peaking background. Below we discuss the sources 
of systematic uncertainties that contribute to our TZ)j!^ 
measurements: 

1. Signal NN shape: in the nominal fit, we use the 
NN PDF from the B signal MC. To estimate the 
related systematics, we refit the data using a sig- 
nal NN PDF extracted from the high purity and 
high statistics B — Dtt RS data, after subtracting 
the residual continuum background contamination 
predicted by the simulation. We set the systematic 



uncertainty to the difference with the nominal fit 
result. 

2. B background NN shape: from a study of generic 
BB MC, it appears that the NN spectra of B back- 
ground events in the niES-^E signal box are similar 
to the signal (but suffer from very low statistics), 
while the NN spectra of background events in an 
enlarged m-ES-^E region differ significantly from 
the signal and show less peaking close to 1. In 
the nominal fit we assumed that both the peaking 
and the non-peaking BB background components 
could be described by the B — > Dtt signal NN 
PDF. To estimate the related systematic error, we 
used BB generic background events selected in a 
AE-niEs enlarged window \AE\ < 200 MeV and 
mEs > 5.20GeV/c2 to build the NN PDF of the 
non-peaking part of the BB background (keeping 
the signal NN PDF to describe the peaking part of 
this background) and repeated the fits, taking the 
difference of the results as the associated systematic 
uncertainty. 

3. Continuum background NN shape: to account for 
possible differences between the simulation and the 
data, we used the NN spectrum from off-peak data 
instead of qq MC (g = u, d, s, c) to model this com- 
ponent. We set the associated systematic uncer- 
tainty to the difference of the two results, but the 
error is dominated by the large statistical uncer- 
tainty on the off-peak data sample. 

4. The shape parameters Cs^^' and Cs*'^'' ^^'^ 
GUS functions describing the suppressed and fa- 
vored BB combinatorial background: in the nomi- 
nal fits, these parameters are fixed to their values as 
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FIG. 3: (color online). Projections on tties (top) and A''A^ (bottom) of the fit results for Dtt (a,d), D'jj^qtt (b,e) and D'^^tv (c,f) 
RS decays, for samples enriched in signal with the requirements NN > 0.94 (rriES projections) or 5.2725 < mES < 5.2875 GeV/c^ 
{NN projections). The points with error bars are data. The curves represent the fit projections for signal plus background 
(solid) and background (dashed). 
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FIG. 4: (color online). Projections on tties (top) and A'^A*' (bottom) of the fit results for Dn (a,d), D'jj^ott (b,e) and D^^^-k (c,f) 
WS decays, for samples enriched in signal with the requirements A*' TV > 0.94 (rriES projections) or 5.2725 < tties < 5.2875 GeV/c^ 
(A'^A'^ projections). The curves represent the fit projections for signal plus background (solid), the sum of all background 
components (dashed), and qq background only (dotted). 



determined from BB simulated events. To account 
for possible disagreement between data and simula- 
tion, we repeated the fits varying these parameters 
in a conservative range. 

5. Peaking component in the B background: we var- 



ied the yield of the peaking component by zblcr, 
where a is either the statistical error from a fit to 
generic BB MC or the uncertainty on the branch- 
ing fraction for known sources of peaking back- 
ground. 
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TABLE IV: Summary of fit results for D^*>-n: 



Mode Dtt ^dtt""" D}y^-K 

Ratio of rates, 7^^^ (10"^) 3.3 ± 0.6 3.2 ± 0.9 2.7 ± 1.4 

Number of signal events Nws 79.8 ± 13.8 28.3 ± 7.7 18.7 ± 9.7 

Number of normalization events Nrs 24662 ± 160 9296 ± 102 7214 ± 105 

B+ ratio of rates, 7^^^+ (10~^) 3.2 ± 0.8 3.5 ± 1.2 4.6 ± 2.2 

B- ratio of rates, 71^^" (10"^) 3.4 ± 0.8 2.9 ± 1.2 1.0 ± 1.8 

Asymmetry A^gj 0.03 ± 0.17 -0.09 ± 0.27 -0.65 ± 0.55 



6. Uncertainty on the number of BB combinatorial 
background events: in the D*tt (and D*K) fits 
where this component has been fixed, we vary it 
by ±25% (the level of agreement between data and 
simulation observed in the Dtt and DK fits) and 
we take the difference with the nominal fit result 
as a systematic uncertainty. 

The resulting systematic uncertainties are listed in Ta- 
ble |Vl We add them in quadrature and quote the results: 

Ud^ = (3.3±0.6±0.4) X 10-^ 
7^^^^o)^ = (3.2 ± 0.9 ± 0.8) X 10-3, 
nia^^^ = (2.7 ± 1.4 ± 2.2) X 10-3, 

where the first uncertainty is statistical and the second 
is systematic. The values of Tlj^^^ are in good agree- 
ment with the world average Rd — fjj — B{D^ — 
K+TT')/B{D" ^ K-n+), Rd = (3.36±0.08) x lO'^ Q. 

A separate fit to B^ and B~ candidates provides a 
measurement of the corresponding asymmetries. We ob- 
tain the following results: 

Aon = 0.03 ±0.17 ±0.04, 
ylp^o)^ - -0.09 ±0.27 ±0.05, 
Alo-y)^ - -0.65 ±0.55 ±0.22, 

where the uncertainties are dominated by the statistical 
error. No significant asymmetry is observed for the 
WS decays. The largest source of systematic uncertainty 
on the D^*^Tr asymmetries is from the uncertainty on the 
B background peaking component. 

B. Results for B D^*^ K 

The results for B — > D^*^K are displayed in Fig. [S] 
(RS modes) and Fig. [5] ( WS modes) . They are summa- 
rized in Table IVII Indications of signals are observed 
in the B DK and in the B -D^^oif WS modes, 
with statistical significances of 2.2a and 2.4cr, respec- 
tively (Fig. [7]) . Accounting for the systematic uncertain- 
ties, the significances become 2.1cr and 2.2a, respectively. 
For B — >■ D*jj K WS, no significant signal is observed. 



TABLE V: Summary of systematic uncertainties on TZ for 
in units of 10"^. 



Source A7^(10-^) A7^(10-^) A7^(10-^) 





Dtt 






Signal NN 


±0.1 


±0.1 


±0.1 


BB background NN 


±0.1 


±0.1 


±0.9 


udsc background NN 


±0.1 


±0.1 


±0.3 


BB comb, bkg shape (mss) 


±0.2 


±0.1 


±0.2 


Peaking background WS 


±0.2 


±0.8 


±2.0 


Peaking background RS 


±0.0 


±0.1 


±0.1 


BB comb, bkg 




±0.0 


±0.4 


Combined 


±0.4 


±0.8 


±2.2 



The systematic uncertainties have been estimated by 
testing different fit models and recomputing TZj^j^, as 
explained in Section IV Al A summary of the different 
systematic uncertainties is given in Table I VIII The un- 
certainties on the NN describing the BB combinato- 
rial background and the uncertainties on the BB peak- 
ing background are the two main contributions. For 
B± ^ DK^, we find for the ratio of the WS to RS 
decay rates 

Udk = (1-1 ± 0.5 ± 0.2) X 10-^ 

Expressed in terms of event yields, the fit result is 
19.4 ± 9.6 ± 3.5 WS events. The resuhs of fits to sep- 
arate B'^ — )■ DK^ and B~ — > DK~ data samples are 
given in Table IVTl Projections of the fits to 5+ and B^ 
data are shown in Figs. [5] and [51 respectively. We fit 
7^+^ = (2.2 ±0.9 ±0.3) X 10-2 for the B+ sample, corre- 
sponding to 19.2 ± 7.9 ± 2.6 events. On the contrary, no 
significant WS signal is observed for the B~ sample, and 
we fit 7^^;^ = (0.2 ± 0.6 ± 0.2) x 10"^. The statistical 
correlation between 7?.^^ and TZ^j^ (or TZdk and Adk) 
is insignificant. 

The systematic errors on the asymmetries are esti- 
mated using the method discussed previously. The main 
systematic error on Adk is from the uncertainty on the 
number of peaking B background events for the WS 
channel. This source contributes lg']^4 to Adk, and 
±0.08 X 10"^ to TZdk, where the changes in the two 
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FIG. 5; (color online). Projections on jties (a, b, c) and A'^A'' (d, e, f) of the fit results for DK (a, d), D'jj^qK (b, e) and 
D}j^K (c, f) RS decays, for samples enriched in signal with the requirements A'^A'' > 0.94 (niES projections) or 5.2725 < mES < 
5.2875 GeV/c^ [NN projections). The points with error bars are data. The curves represent the fit projections for signal plus 
background (solid) and background (dashed). 




FIG. 6: (color online). Projections on jties (a, b, c) and A'^TV (d, e, f) of the fit results for DK (a, d), D^^qK (d, e) and D'jj^K 
(c, f) WS decays, for samples enriched in signal with the requirements A'"A'^ > 0.94 (tties projections) or 5.2725 < mES < 
5.2875 GeV/c^ (A^A'^ projections). The points with error bars are data. The curves represent the fit projections for signal plus 
background (solid), the sum of all background components (dashed), and qq background only (dotted). 



quantities are 100% negatively correlated (increasing the 
peaking background increases Adk but decreases TZdk)- 
The other sources of systematic uncertainty considered in 
Table IVlIl are 100% correlated between 72.+ and 7?.~, and 
mostly cancel in the asymmetry calculation. By com- 



paring the number of and B~ events reconstructed 
in the [i^+7r+]£)7r+ analysis, where no significant asym- 
metry is expected, the uncertainty due to the detector 
charge asymmetry is estimated to be below the 1% level. 
Finally, we also account for a possible asymmetry of the 
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TABLE VI: Summary of fit results for D'^*^K 



Mode DJK Dp^oK D*d^,K 

Ratio of rates, 7^^]^ (10"^) 11.1 ± 5.5 17.6 ± 9.3 13 ± 14 

No. of signal events Nws 19.4 ± 9.6 10.3 ± 5.5 5.9 ± 6.4 

No. of normalization events Nrs 1755 ± 48 587 ± 28 455 ± 29 

B+ Ratio of rates, 7^^]^ (10~^) 21.9 ± 9.0 4.9 ±7.9 9 ± 16 

B- Ratio of rates, 71^]^" (10"^) 1.7 ± 5.9 37 ±18 19 ± 23 

Asymmetry y^^]^ -0.86 ± 0.47 0.77 ± 0.35 0.36 ± 0.94 




0.01 0.02 0.03 0.04 0.05 0.01 0.02 0.03 0.04 0.05 0.01 0.02 0.03 0.04 0.05 



FIG. 7: Negative log-likelihood variation vs 7^^]^ for -> DK^ (left), B^ D*jj-„oK^ (center) and B^ D^^K^ (right). 
Systematic uncertainties are not included. 



charmless — K^K^ir^ peaking background. The 
asymmetry of this background has been measured to be 
± 10% [13] and we estimate the corresponding system- 
atic uncertainty by assuming a ±10% asymmetry of this 
background. The final result for the asymmetry is: 

=-0.86 ±0.47 tlfa ■ 



TABLE VII: Summary of systematic uncertainties on TZ for 
L>(*)A', in units of 10"^ 



Error source S^^TF"^y^^^TO^'^y^^^T(r^ 





DK 






Signal NN 


±0.1 


±0.1 


±0.3 


BB background NN 


±0.1 


±0.3 


±0.4 


qq background A'^A'' 


±0.1 


±0.1 


±0.1 


BB comb, bkg shape (tties) 


±0.1 


±0.1 


±0.1 


Peaking background WS 


±0.2 


±0.3 


±0.6 


Peaking background RS 


±0.0 


±0.1 


±0.1 


Floating BB comb, bkg 




±0.1 


±0.2 


Combined 


±0.2 


±0.4 


±0.8 



For D*jj^oK^, we find for the ratio of the WS 

to RS decay rates 

T^D^^')K = (1-8 ± 0.9 ± 0.4) X 10-2. 



Expressed in terms of event yields, the fit result is 
10.3 ± 5.5 ± 2.4 WS events. The resuhs of fits to sep- 
arate _B+ — ^ D*K+ and B^ — > D*K^ data samples are 
given in Table IVTl Projections of the fits to 5+ and B~ 
data are shown in Figs. [S] and El respectively. Wc find 
'^{DnO)K = (3-7 ± 1.8 ± 0.9) X 10-2 for the B- sample, 
corresponding to 10.2 ±4.8 ±2.4 events. On the contrary, 
no significant WS signal is observed for the i?+ sample, 
and we find 7^*+^o)^^. = (0.5 ±0.8 ±0.3) x lO^^. The sys- 
tematic errors are estimated using the same method as 
for B^ — > DK^, separately for and B~ events. The 
main systematic error on the asymmetry -^(^ijrO)^ from 
the uncertainty on the number of peaking B background 
events for the WS channel. This source contributes ±0.09 
to A*^^o^K, and tO.3 x 10"^ to 7^p^o)^, where the 
two quantities are anti-correlated. The other sources of 
systematic uncertainties mostly cancel in the asymmetry 
calculation, because they induce relative changes on 7?.*+ 
and TZ*~ which are 100% correlated. The final result for 
the asymmetry is: 

Ala-^a)K = +0.77 ± 0.35 ± 0.12. 

The asymmetry for D'^^qK has the opposite sign to the 
asymmetry for DK, in agreement with the shift of ap- 
proximately 180° between Sb and S'^ suggested by the 
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measurements of Refs. [EQ- 

For B — > D*jj^K , we have no significant signal and fit 

= (1-3 ± 1-4 ± 0.8) X 10-2. 

Expressed in terms of event yields, this result corresponds 
to 5.9 ± 6.4 ± 3.2 events D^^^K WS. We fit 211 ± 19 RS 
B- events and 244 ± 20 RS B+ events, and find for the 
WS to RS ratios = (1.9 ± 2.3 ± 1.2) x IQ-^ and 

^(i)7)i^ = (0-9 ± 1-6 ± 0.7) X 10-2. The corresponding 
asymmetry is 



TABLE VIII: Constraints on from the combined B — > 

[Kii\j^{,)K~ ADS measurements. 



Parameter 


1 (7 meas. 


90% C.L. upper limit 




(9.5+li)% 


< 16.7% 


r*g from 








(13.11^:?)% 


< 19.5% 




(12.0tl«;°)% 


< 24.5% 


all D*" decays 


(9.61^;?)% 


< 15.0% 



^(D7)K = +0-36±0.94l[5;25. 



VI. DISCUSSION 



We use the B D^*'>K analysis results and a fre- 
quentist statistical approach [Tsj to extract information 

on and r^K In this technique a is calculated us- 
ing the differences between the measured and theoretical 
values (including systematic errors) of the various ADS 
quantities from Eqs. (P), ^ and ([5]). We assume Gaus- 
sian measurement uncertainties. This assumption was 
checked to be valid and conservative at low tb values 
with a fuU frequentist approach i5|. For B~ — >• DK~, we 
have for instance 



X = 









+ (T^DK- 




.rD,5D)f /cf\- 


+ ir^D^ - 






+ (4"' - 




(9) 



summarized in Tab. IVIIII For B^ — > [Ktt\i:,K^ , we find 
the minimum at = (9.5^4 it)%. This leads to the 
upper limit: rs < 16.7% at 90% C.L., to be compared 
to < 23% at 90% C.L. for the previous ADS analysis 
as performed by BaBAR [i| with 232 x 10^ BB pairs, and 
to rs < 19% at 90% C.L. for the corresponding ADS^ 
analysis as performed by Belle ^ with 657 x 10^ BB 
pairs. We exclude rs — with a C.L. of 95.3%. Similarly, 
for B- [Kit]d-K- we find r*g = (9.6^^;^)%. This 
leads to the upper limit: < 15.0% at 90% C.L., to 
be compared to < 16% at 90% C.L. for the previous 
BABAR ADS analysis [l]. We exclude = with a C.L. 
of 83.9%. 

Using the above procedure we also determine the 2D 
confidence intervals for 7 vs (5^' shown in Figs. 1111 and 
[T^ Choosing the solution with < 7 < 180° favors a 
positive sign for the strong phase 5b [Adk < 0), and 
a negative sign for the strong phase Sg {■^*dttO)k > '^)- 
This result is in good agreement with the values of the 
strong phases determined in Refs. 0, 01- Finally, FigfT51 
shows the C.L. curve as a function of 7 when combining 
the DK and D*K results. 



where Tl^^\rB,l,5B,rD,5D) is given by Eq. (P), and 
where the two last terms constrain and 5d to the 
values r^-* and J^-* of Ref. [s"] within their errors ct^ 
and as. The choice of (7^^^, T^dk) J^ather than {TZdk, 
Adk) is motivated by the fact that the set of vari- 
ables (TZdk, Adk) is not well-behaved (the uncertainty 
on Adk depends on the central value of TZdk), while 
i'^'oK' "^dk) ^""^ statistically independent observ- 
ables. In the same way, the two pairs of ADS observables 
i'^*inO)K^ T^lDnO)K) and (7^(^.^)^, 7^(^^)^) are used to 
extract r^, while accounting for the relative phase differ- 
ence in the two D* decays @. We allow < < 1, 
-180° < 7 < 180°, and -180° < < 180°. The mini- 
mum of the for the rg \ 7, dg \ td, and 5d parameter 
space is calculated first (x^in)- We then scan the range 



of minimizing the (Xm) by varying 5^'', 7, td, 
and 5d- A confidence level (C.L.) for tb is calculated 



(*) 



using Ax 



and one degree of freedom. 



The results of this procedure are shown in Fig. [TUlfor 
the C.L. curve as a function of . The results are 



VII. SUMMARY 

In summary, using a data sample of 467 million BB 
pairs, we present an updated search of the decays B~ — 
D^*'>K~ where the neutral D meson decays into the 
K^TT^ final state (WS). The analysis method is first ap- 
plied to B~ — D^*'>'K~, where the D decays into the 
Cabibbo favored (K^ir^) and doubly suppressed modes 
(K+TT-). We measure TZd^ = (3.3 ± 0.6 ± 0.4) x 10-^, 
T^hrro)^ = (3-2 ± 0.9 ± 0.8) X 10-3 and 7^^^^), = 
(2.7 ± 1.4 ± 2.2) X 10-3, in good agreement with the ra- 
tio Rd of the suppressed to favored D° — > Kn decay 
rates, Rd = (3.36 ± 0.08) x lO'^ Q. Both the branch- 
ing fraction ratios and the CP asymmetries measured 
for those modes, Adtt = (3 ± 17 ± 4) x 10-^, A*^^^o-)^ = 
(9±27±5) X 10-2 and ^(^^)^ = (65±55 tH) x IO-2, are 
consistent with the expectations discussed in Section H] 

We see indications of signals for the B — >■ DK and 
B — >■ D*jj^oK wrong-sign modes, with significances of 
2.1(7 and 2.2(7, respectively. The ratios of the WS to RS 
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FIG. 8: (color online). Projections on rriES (a, b, c) and A^A'^ (d, e, f) of the fit results for DK^ (a, d), D'^^^oK^ (b, 
e) and D*jj^K^ (c, f) WS decays, for samples enriched in signal with the requirements A''A'^ > 0.94 (mES projections) or 
5.2725 < rriES < 5.2875 GeV/c^ (A'^A'^ projections). The points with error bars are data. The curves represent the fit projections 
for signal plus background (solid), the sum of all background components (dashed), and qq background only (dotted). 
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FIG. 9: (color online). Projections on tties (a, b, c) and TVA'^ (d, e, f) of the fit results for DK (a, d), D^^qK (b, 
e) and D^^^K' (c, f) WS decays, for samples enriched in signal with the requirements A''A' > 0.94 (mES projections) or 
5.2725 < rriES < 5.2875 GeV/c^ {NN projections). The points with error bars are data. The curves represent the fit projections 
for signal plus background (solid), the sum of all background components (dashed), and qq background only (dotted). 



branching fractions are measured to be TZdk = (1-1 ± 
0.5±0.2) X 10-2 and 7^^£,^o)K = (1.8±0.9±0.4) x 10"^ for 
B DK and B — > D*jj^oK, respectively. The separate 
measurements of 'R.jjj^ for B^ and B" events indicates 
large CP asymmetries, with Adk = —0.86 ± 0.47 ^o'}g 



for B DK and -4(^^o)^ = +0.77 ± 0.35 ± 0.12 for 
B ^ D*K, D* ^ Dt:°. For the B D*jj^K WS mode, 
we see no statistically significant evidence of a signal. We 
measure 'R-(D.y)K = (l-3± 1.4±0.8) x 10"^ and A^jj^^j^ = 
+0.36 ± 0.94 toll. These results are used to extract the 
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FIG. 10: (color online). Constraints on from the com- 
bined — > [K 7v]jj(=i.) K ADS measurements. The solid 
(dotted) curve shows the 1 minus the confidence level to ex- 



clude the abscissa value as a function of r 



The horizontal 



lines show the exclusion limits at the 1 and 2 standard devi- 
ation levels. 
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FIG. 12: (color online). One minus confidence level isocon- 
tours on 7 vs 5% from the combined — >■ [Kn]D* ADS 
measurements. 



1-CL 





150 




100 




50 











0) 




2, 















-50 




-100 




-150 



1 






1 1 1 1 1 1 1 1 1 1 1 1 1 1 





-150 -100 



50 100 150 



So (deg) 



1.0 



0.8 



0.7 



0.6 



0.5 



0.4 



0.3 



0.2 



0.1 



0.0 



FIG. 11: (color online). One minus confidence level isocon- 
tours on 7 vs 5b from the B~ — >■ [Kn]DK~ ADS measure- 
ment. 













2 a 




- 1 . . , , 1 , , , , 1 , , , . 1 , . . , 1 . , , , 1 , , - 


-150 


-100 


-50 





50 


100 


150 



Y (deg) 



FIG. 13: (color online). Constraints on 7 from the combined 
B~ — )■ D^*^[K^TT~]K~ ADS measurements. The solid curve 
shows the (1-C.L.) to exclude the abscissa value. The hori- 
zontal lines show the exclusion limits at the 1 and 2 standard 
deviation levels. 



following constraints on Tg . 

TB = (9.5±5;1)%, 

r*B = (9.6l^5j)%- 

Assuming < 7 < 180°, we also extract constraints 
on the strong phases 5^ , in good agreement with other 
measurements Ref. [1, Q . 
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